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Abstract 

New Brillouin scattering measurements of velocity and attenuation of sound in the hypersonic 
regime are presented. The data are analyzed together with the literature results at sonic and 
ultrasonic frequencies. As usual, thermally activated relaxation of structural entities describes the 
attenuation at sonic and ultrasonic frequencies. As already shown in vitreous silica, we conclude 
that the damping by network viscosity, resulting from relaxation of thermal phonons, must be 
taken into account to describe the attenuation at hypersonic frequencies. In addition, the bare 
velocity obtained by subtracting to the experimental data the effect of the two above mechanisms 
is constant for temperatures below 250 K, but increases almost linearly above, up to the glass 
transition temperature. This might indicate the presence of a progressive local polyamorphic 
transition, as already suggested for vitreous silica. 



I. INTRODUCTION 



In insulating solids, high frequency atomic vibrations are responsible for heat trans- 
port and thus determine the thermal properties. While the mechanisms leading to phonon 
damping in crystals are rather well understood, those responsible for sound attenuation in 
non- crystalline solids are still a matter of large interest. The main problem is that, due to 
the lack of periodicity, the exact nature of local and of collective vibrations is unknown. 

At very low temperatures, it has been demonstrated that the dominant damping mech- 
anism is a coupling of acoustic vibrations with tunneling units 1 '-. While the microscopic 
nature of those systems is not known, the phenomenological model involving resonance and 
relaxation of two-level tunneling systems (TLS) explains the temperature (T) and frequency 
(Q) dependence of sound velocity v and attenuation a2&. In a large T domain above 10 K 
the sound attenuation is dominated by a large broad peak which shifts to high T with 
increasing frequency. This peak has been assigned to the coupling of acoustic vibrations 
with structural entities described as defects moving in asymmetric double-well potentials by 
thermally activated relaxation (TAR)£~— . This explains the T-dependence of a observed in 
most sonic and ultrasonic experiments. However, the T-dependence of the sound velocity in 
glasses such as AS2S3 and As2Ses 8 cannot be explained by this model: the behavior observed 
there is very similar to the decrease of sound velocity in crystals, where it is assigned to 
relaxation of thermal phonons by anharmonic interactions*^. This suggests that, even if 
the nature of thermal modes in glasses is unknown, and certainly different from propagating 
plane- wave phonons, a similar anharmonic mechanism can be presentii^. In fact, we have 
observed that Brillouin scattering measurements of a in vitreous silica at frequencies in the 
hypersonic range cannot be described satisfactorily by the TAR mechanism alone^. The 
excess to the TAR contribution was assigned to anharmonicity. The measurements were 
extended to frequencies up to 400 GHz^. The results are well described by the sum of the 
contributions of the three above mechanisms, TLS, TAR and anharmonicity. Due to the 
different T and Q dependence of those mechanisms, anharmonicity tends to dominate above 
100 GHz. In vitreous silica, t>-Si02, we have performed a simultaneous fit of measurements 
of a and v at low T in a large frequency domain, which allowed to extract the parameters 
of the models. Calculating the velocity dispersion up to high T and subtracting the con- 
tributions given by TAR and anharmonicity from the experimental results, one expects to 
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get a constant bare velocity. In fact, we observed that the bare velocity shows an almost 
linear increase above 100 K. This was explained by a temperature dependent microscopic 
structural change leading to T-hardening^"— . 

A well-known characteristic of glasses is the plateau observed in the thermal conductiv- 
ity in the 5-10 K range 18 . This plateau suggests that a strong damping mechanism limits 
the mean free path of the acoustic waves responsible for thermal transfer—. In this tem- 
perature range, the frequency of dominating thermal acoustic waves is around 1 THz 20 . A 
mechanism leading to this strong damping must be invoked in addition to those discussed 
above. Using Brillouin scattering of X-rays to investigate acoustic waves in the THz range, 
we have given evidence for an attenuation increasing as Q 4 in two glasses leading to a Ioffe- 
Regel crossover—. Above the frequency of this crossover, plane acoustic waves can no longer 
propagate^. 

The association of the mechanisms discussed above gives a satisfactory description of 
most of the a and v measurements in vitreous silica. The purpose of the present paper 
is to extend this description to vitreous germanium oxide (f-GeO^), another tetrahedrally 
coordinated glass. 

II. MATERIAL AND METHODS 

The samples used were kindly provided by Dr. J. C. Lasjaunias, from the Centre de 
Recherche sur les Tres Basses Temperatures in Grenoble, France. Their optical transparency 
and homogeneity are excellent. A rough estimate by infrared absorption at 3560 cm -1 gives 
a value of ~ 100 ppm for the OH" content. We have observed the influence of the thermal 
history on a measurements in the hypersonic range: differences of about 20% were observed 
between fast quenched and annealed samples. The latter have been used for this study. 
The velocity and attenuation of acoustic waves in the hypersonic frequency range were 
measured by Brillouin scattering (BS). The spectrometer used for the experiments has been 
described elsewhere^. Its resolving power of about 2xl0 7 is large enough to allow precise 
measurements of Brillouin frequency shifts and linewidths. The 514.5 nm line of a single- 
frequency argon-ion laser was used. The measurements were performed in the backscattering 
geometry were only longitudinal modes are active in BS for isotropic media. 

In order to reach the high resolving power needed with a sufficient luminosity, we use a 
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FIG. 1. Typical spectra registered at three temperatures. A full interference order of the confocal 
FP analyzer is shown. The Stokes Brillouin peaks at interference order p are indicated by the 
arrows whereas dashed lines indicate the elastic peaks at orders p+15 and p+16. The red lines 
result from fits explained in the text. 

confocal Fabry-Perot (FP) interferometer with a free spectral range of 1.497 GHz. Some 
typical spectra are shown in Fig. [TJ The Brillouin frequency shift is about 24 GHz, which 
means that the elastic and the Brillouin lines are observed at different interference orders. 
The frequency-shift scale on the figure is that for the Brillouin line and the arrows indicate 
the central frequency of the latter. The excellent signal-to-noise ratio observed for the spectra 
is reached in about 20 minutes at 5.5 K whereas 5 minutes are sufficient at room-T. From 
this figure, it is clear that both the frequency shift and the linewidth can be measured with 
a high accuracy in the whole temperature range. To extract these parameters, the Brillouin 
spectra are fitted to a damped harmonic oscillator convoluted with the instrumental profile 
taken from the elastic line. The small parasitic broadening due to the aperture of collection 
(~ 40 mrad) is taken into account. 

III. RESULTS AND DATA ANALYSIS 

The results for T from 5.5 K to 900 K are shown in Fig. [2j From the lowest T, the 
frequency shift SQb/^ decreases strongly to a minimum around 300 K and increases above, 
up to the glass transition temperature T g ~ 800 K. At the latter T, a cusp is observed, 
as usual in glasses. The fullwidth at half maximum of the Brillouin line, r/27r, which is 
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FIG. 2. Brillouin frequency shift SVLb/^tt (squares, left scale) and fullwidth T/2ir (triangles, right 
scale) as a function of temperature. The vertical dashed line marks the glass transition temperature 
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very small at low T, increases rapidly with increasing T up to a maximum around 300 K, 
decreases above up to ~ 700 K, and finally shows a plateau at large T, without noticeable 
change at T g . 

From the values of 5Qb, the velocity v can be calculated using 

5Qb 2nv 



2tt 



Ao 



sin 0/2, 



(1) 



where n is the refractive index, Ao the wavelength of the incident light in vacuum and 9 
the scattering angle. The values of n = 1.610 at 300 K and its T-dependence were taken 
from the literature^. The attenuation a is obtained from T measurements, a = Y/v. For 
comparison of attenuation measurements in a large Q domain it is more convenient to use 
the internal friction: 



Q- 1 



5Q 



B 



(2) 



A previous Brillouin scattering measurement of v and Q _1 in v-Ge02 was performed 
earlier by Hertling et al.— in the T-range from 50 K to 750 K. Our results which are much 
more accurate agree with these measurements. 



IV. DISCUSSION 



In order to elucidate the contribution of the different damping mechanisms in vitreous 
germania, it is necessary to analyze simultaneously data measured in large T and Q domains. 
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FIG. 3. Internal friction of u-GeO^ as a function of temperature at several frequencies: 6.3 kHz», 
50 kHz 26 , 7.8 MHz 2 ! 10 MHz^, 12 MHz^E, 20 MHz^, 21.2 MHz 2 !, and 50 MHz^. 

In comparison to the situation in vitreous silica, the data in t>-GeO"2 are scarce. The results 
of sonic and ultrasonic measurements collected from the literature are shown in Fig. [31 As 
in most other glasses, a strong broad peak is observed in the Q -1 measurements over the T, 
Q domain shown. This indicates that the TAR mechanism is the main origin of the sound 
damping there. From our study of vitreous silica and also from preliminary calculations of 
the possible anharmonic damping in t>-GeO"2, we know that this contribution is negligible in 
this T, Q domain. This is also the case for Q" 1 from the lower T up to 50 K in Brillouin 
scattering experiments. The temperature shift of the maximum with increasing Q observed 
in Fig. |3] allows to evaluate an average activation energy Vq ~ 2400 K, as already shown by 
Hertling et al.—. It is more difficult to compare the amplitude of the curves, as measurements 
at almost the same Q give appreciably different values. In particular, the measurements at 
7.8, 12 and 21.2 MHz have almost the same amplitude, but are lower by a factor of about 
1.7 than those at 10, 20 and 50 MHz. The comparison of amplitude from measurements 
performed on glass samples of different origin is generally difficult, as the thermal history 
and OH content have an appreciable influence on the results. On the other hand, absolute 
calibration of Q~ x measurements in ultrasonic experiments is always a difficult task. It is 
thus impossible to fix in the fit the same amplitude coefficient for the TAR contribution 
from the experiments shown in Fig. |3j 

The thermally activated defects are described by a distribution P(A, V) of asymmetric 
double-well potentials with a barrier height V and an asymmetry A. The contribution of 
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FIG. 4. Internal friction of u-GeG^ in the sonic and ultrasonic frequency ranges fitted to the sum 
of TAR and TLS contributions as discussed in the text. Thin black lines show the small TLS 
contributions extracted from Refs. 25 and |28|. 



the TAR mechanism to Q 1 (see e.g. Ref. 13JJ) , Q T a R , is given by 

2 poo PCO A Q 

where 7 is a deformation potential, and p is the mass density. In this equation, r is the 
relaxation time for hopping within the double well. It is given by 

r = r exp^sech^, (4) 

where Tq is the inverse of an attempt frequency^. 

We have made a common fit of all Q~ l measurements shown in Fig. 131 together with our 
Brillouin data at T below 50 K. As a consequence of the above discussion on the amplitudes, 
we have fixed a common amplitude coefficient for measurements at 10, 20, 50 MHz and 
24 GHz. This coefficient was left free for measurements at the other frequencies. 




200 400 600 800 

Temperature (K) 



FIG. 5. (triangles) Internal friction of v-Ge02 at hypersonic frequencies derived from the Brillouin 
linewidths. The dashed line shows the TAR contribution calculated from Eq. ([3]) whereas the thin 
solid line is the TLS contribution. The circles represent the internal friction observed in Brillouin 
scattering after subtraction of Q^AR anc ^ ^tls - ^ ne dotted line is the calculated anharmonic 
contribution Q^nh as ex pl arne d in the text. The solid line is the sum of these three contributions. 

For the distribution P(A, V), our first trial was a gaussian distribution for both param- 
eters as in Ref. [l3|. The result of the fits disagrees with the experiments, giving a profile 
which is much broader than that observed. A better description was obtained by taking a 
sharper cut-off at Vq, i.e. an analytic form proportional to exp(— (V 2 /2Vq) 2 ) for the distri- 
bution of barriers. Furthermore, the Q~ l calculation at the frequencies of the experiments 
shown in Fig. [3] is hardly sensitive to the cut-off Ao in the distribution of asymmetries as the 
corresponding relaxation times are too short to influence the acoustic damping. In contrast, 
the velocity change 5v , given by 



5v 
v 



1 7 2 



dA / dVP{A, V) sech 2 



A 



1 



-2 : 



(5) 
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is sensitive to the integral over all defects having short relaxation times. For this reason, we 
have also included the values of 5v measured in the temperature region where the anharmonic 
contribution can be neglected, i.e. below 100 K. The results of the fits in the sonic and 
ultrasonic regimes are shown in Fig. HJ One finds Vq ~ 2460 K, To ~ 0.6 xlO -14 s and Ao — 
180 K. The two first values are in good agreement with those found in Ref. bj. The result 
for A cannot be compared with Ref. [2J as the authors used a distribution with a constant 



density. Considering the differences in profiles between experiments taken at almost same 
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FIG. 6. (a) Qtar sca l e d by A(T)/^ h . The solid line is obtained with Eq. ([6]) using the fitted 
values of Tth{T) in (b). (b) Temperature dependence of the thermal time r-th- The solid line is a 
T~ 2 power law adjusted to the data. 

frequencies, the agreement is satisfactory. It must be noted that the small contribution 
of TLS to Q^ 1 was calculated from parameters extracted from the literature 2 ^ 1 ^. This 
contribution is shown in Fig. H] as a thin black line and was added to the calculated TAR 
value to allow comparison with experiment. 

From the parameters determined for the TAR mechanism, it is now possible to calculate 
the contribution of this process to Q _1 at Brillouin scattering frequencies. The result is 
shown as a dashed line in Fig. [5j It is clearly much below the measured values. By subtract- 
ing Qtar an d Qtls fr° m the measurement, one gets the anharmonic contribution, Q^nhj 
shown by circles in Fig. [5j One finds that Qanh increases with increasing T up to about 
300 K and is nearly constant above up to T g . 

The anharmonic contribution in crystals was calculated by Maris—. To analyze Qanh h 1 
glasses, we follow a procedure described and discussed in detail in Ref. 

Qanh i s gi ven 



by: 

Qanh = A- ■ 02 ^ 2 , (6) 



1 + Wrl 

with 

2 C v Tv 

A = 7th it— r- ( 7 ) 



'v 

2pvl 



Here, C v is the specific heat per unit volume and vd is the Debye velocity, whereas r t h and 
7 t 2 h are the average lifetime and the mean-squared average Griineisen parameter of thermal 
phonons, respectively. Except for the value of 7 2 h , all the quantities in A are known for v- 
Ge0 2 according to Eq. (jTJ). The low-T values of C v were taken from Ref. [l8|, complemented 




FIG. 7. (squares) Velocity changes of u-GeC^ in the sonic and ultrasonic frequency ranges. The 
solid lines are the sum of TAR (dashed) and anharmonic (dotted) contributions as discussed in the 
text. 



by the high-T values of Ref. |33j- We used vd ~ 2800 m/s and p ~ 3.65 g/cm 3 . In Fig. [6^, 
we have plotted Tth^ANHM as a function of T. From Eq. ([6]), one expects a maximum 
of this quantity at a value 7th/2, when f2r t h = 1. From the points in Fig. EK, we observe 
that Tth^5 anh / increases when T decreases up to a shallow maximum around 0.4. From 
this value, we extract an evaluation of 7 t 2 h ~ 0.8. It is now possible to calculate 7~th &s a 
function of T. The result is shown in Fig. [6]d. Except for the values below 200 K for which 
the uncertainty is large, the results agree with a variation r t h oc T -2 . The T dependence 
of T^QanhM can now be calculated: the result is shown as a solid line in Fig. [6k With 
the same parameters, it is possible to calculate Qanh fr° m Eq. ([6]). The dotted line in 
Fig. [5] shows the result, which is in good agreement with the estimation of this contribution 
extracted from the experiment. Finally, <3tls> Qtar an d Qanh are added: the solid line in 
Fig. [5] is obtained, in excellent agreement with the experiment. 

From the known parameters for TAR and anharmonicity, it is possible to calculate the 
velocity changes produced by these mechanisms, <5t>TAR an d <5^anh, respectively (see Eq. ([5]) 



for TAR and Eq. (16) of Ref. [l_3j for anharmonicity). The results are shown in Fig. [7] where 
they are compared to the experimental values at sonic and ultrasonic frequencies. It must 
be noted that 5v anh is nearly constant below 100 K, and that Svtar shows a deep minimum 
around 200 K, while 5v anh decreases linearly in this region. On the three figures [T^-c, we 
remark that the slope of 5v(T) below the minimum is well reproduced. This slope is directly 
related to A , the upper cut-off in the distribution of asymmetries. This indicates that 
only relaxing entities with asymmetry below about 180 K are present in f-Ge02- On the 
other hand, the shallow minimum observed in the measurement at 10 MHz (Fig. [7]d) is well 
reproduced. In contrast, the calculation disagrees with the experiment at 20 MHz. It must 
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FIG. 8. (squares) Velocity changes in f-GeO"2 at Brillouin scattering frequencies. The dashed 
(dotted) line shows the calculated TAR (anharmonic) contribution whereas the solid line is the 
sum of these two mechanisms. The dashed-dotted line represents the relative changes of the 
calculated bare velocity Voo. 

be noted that only a very small change of 5v is expected when the frequency varies from 
10 MHz to 20 MHz as shown by the solid lines in Fig. [7)>c. Thus, we suggest that the 
disagreement at 20 MHz originates from uncertainties in the measurement. 

Fig. [8] shows the experimental and calculated values for 5v/v at Brillouin scattering 
frequencies, around 24 GHz. The velocity minimum observed here is much less deep by a 
factor of ~ 2 than at ultrasonic frequencies. This is well described by the calculation, which 
is in excellent agreement with the experiment from the lowest T to ~ 300 K. From this 
observation, and from the above remarks on the velocity slope, we conclude that the density 
of relaxing defects used for the fit is satisfactory. Finally, subtracting 5v tar and Svanh from 
the experimental values, the bare velocity is obtained. In principle, one expects that 
is independent from T. We have plotted {5v/v) ao versus T as a dashed-dotted line in Fig. 
While this quantity is constant from to ~ 250 K, it shows a nearly linear increase above, 
up to T g . A similar behavior was already observed in f-SiOV 3 . 

To summarize, we have shown that a satisfactory description of velocity and attenuation 
measurements in a large T, Q domain can be obtain by taking into account both TAR and 
anharmonicity mechanisms. We have shown that a rather low cut-off in the asymmetry 
distribution of TAR is needed for a good description of the velocity slope at low T. While 
the effect of anharmonicity on Q~ x is negligible at ultrasonic frequencies, the contributions 
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of both processes at Brillouin scattering frequencies have almost the same amplitude. This 
might explain why the velocity measurements in Re f. I24J and the Q~ l values at higher fre- 
quencies measured in UV Brillouin scattering in Ref. Qcould not be described satisfactorily 
on the basis of the TAR model alone. 



V. CONCLUSIONS 

We have shown that the description of acoustic velocity and attenuation already estab- 
lished in f-Si02 is also appropriate for v-Ge0 2 - In both glasses, the effect of anharmonicity 
on the attenuation is negligible in the sonic and ultrasonic frequency range. The changes 
with frequency of the peak in Q~ l observed in this region gives a value Vq — 2460 K for the 
activation energy in t>-Ge02, while Vq ~ 660 K was found in t>-Si02, in agreement with the 
fact that the free-volume is smaller in t>-Ge02 than in f-SiOv The anharmonic contribution 
to Q^ 1 becomes similar to the effect of TAR in the hypersonic regime. The effect of anhar- 
monicity is negligible on the variation with T of the sound velocity below 100 K. The latter 
can be described by assuming a cut-off in the asymmetry of TAR defects. This cut-off was 
found equal to ~ 90 K in t>-SiO"2, while we find Ao — 180 K in v-GeO^. 

In f-Si02, it was shown that, when the velocity changes induced by TAR and anhar- 
monicity are taken into account, the bare velocity is constant at low T, then increases 
linearly above ~ 150 K. When this model is applied to t>-Ge02, we find the same behavior: 
Voo is constant up to 200 K and increases linearly above 300 K. In t>-Si02, this anomalous 
temperature hardening was assigned to changes in the ring conformation induced by flips 
of Si— O— Si bonds^. Our results indicate that a similar mechanism might be present in 
w-Ge02- This could be an important matter to be studied by atomic simulations. 
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